Fluorescence intermittency in blinking quantum dots: renewal or slow modulation? 
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Oh We study time series produced by the blinking quantum dots, by means of an aging experiment, 

and we examine the results of this experiment in the light of two distinct approaches to complexity, 
renewal and slow modulation. We find that the renewal approach fits the result of the aging exper- 

£f>l . iment, while the slow modulation perspective does not. We make also an attempt at establishing 

^s^j ' the existence of an intermediate condition. 
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I. INTRODUCTION 



^ ■ In the last few years there has been an increasing interest for the intermittent fluorescence of new nanomaterials 
0,01, which for this reason have been called Q Blinking Quantum Dots (BQD). The theoretical interest for these 
new materials is due, to a great extent, to the properties of their power spectrum. Pclton, Grier and Guyot-Sionnest 
0, as well as Chung and Bawendi have studied macroscopic samples of quantum dots, and have proved that they 
generate 1// noise. Thus, understanding these new materials is a problem directly connected with the origin of the 
1// noise, an issue that according to the proponents of Self Organized Criticality (SOC) should be settled using their 
theoretical perspective The search for the dynamical origin of complexity, meant to be a statistical condition 
1 departing from the Poisson prescription, is a hot topic, and recently another proposal of approach to complexity, 
'— called super statistics, was made by Beck and Cohen pfl (see also Ref.Q)- On the other hand, there is an increasing 
conviction 0, 0, 0, 0] that anomalous diffusion, and especially sub-diffusion, might find a satisfactory theoretical 
framework in the subordination perspective. In this paper we shall illustrate the subordination perspective with 
qq 1 physical arguments, and shall refer to this theoretical approach simply as subordination. Here we limit ourselves to 
saying that according to the original work of Montroll and Weiss |13). the time distance between two consecutive jump 
^sO , events, generating the walker diffusion, is not fixed but is randomly drawn from a distribution of waiting times. This 
ON ' is a case used to generate a new form of diffusion (anomalous), through subordination to the ordinary diffusion. The 
subordination perspective is also be applied to generate anomalous fluctuation-dissipation processes. In this case the 
original events are processes where diffusion is balanced by friction. By assuming that the time distance between two 
consecutive events of this kind is not fixed, we derive a non-ordinary fluctuation-dissipation process subordinated to 
the ordinary one [l4| . 

Subordination yielding anomalous transport can be imagined as a theoretical perspective where the collisional 
events generating randomness are rare and separated by large time intervals with a non-exponential distribution. This 
important property requires, in turn, a physical explanation that hopefully might relate the deviation from exponential 
distribution to the cooperation among the elementary constituents of the system. The subordination perspective is 
based on the assumption that the phenomenon under study is a renewal one, insofar the jump occurrence is assumed 
^ I to reset to zero the system's memory. 

^ . In the case of BQD fluorescence, many researchers have been studying the microscopic origin of intermittency, and, 
I^j ' on top of that, of why the intermittency of these materials is characterized by non-Poisson distribution of times of 
sojourn in the "on", with the system emitting light, and "off" state, with the system emitting no light p^llTfllr^. IT^ | . 
To first sight, all these papers seem to share the same theoretical perspective as superstatistics. According to Cohen 
the process yielding the deviation from the canonical distribution is realized by the system as follows. The particle 
resides in a cell in local equilibrium for a time extended enough as to reach the corresponding canonical equilibrium, 



'Electronic address: sbianco@unt.edu 
t Electronic address: grigo@unt.edu 



tElectronic address: paradisi@le.isac.cnr.it 



2 



then it moves to another condition of local equilibrium, and so on, in such a way that the resulting distribution will 
turn out to be distinctly non-exponential, insofar as the superposition of many exponential functions with different 
probabilistic weights is not an exponential function. The sojourn in a condition of local equilibrium must be extended 
in time enough for the system to adapt itself to this condition. For this reason superstatistics is equivalent to a form 
of very slow modulation. 

The model discussed in Ref. 0] might be thought of as a form of superstatistics. In fact, these authors propose a 
two-state model with barrier height or barrier width fluctuation. This automatically yields a waiting time distribution 
as a sum of infinitely many exponential functions, one exponential for each of the barrier heights or widths. However, 
as the authors point out, these fluctuations cannot be slow. In fact, the case of slow fluctuations would produce a 
correlation between two successive on-times or two successive off-times, while the analysis of Ref. 0] shows that no 
significant correlation of this kind exists. Similar remarks apply to the model of Refs. uji W$ • I n a U these models 
there are traps with exponential waiting times distribution densities. However, the system resides in each exponential 
trap corresponding to a given state, for instance the state "off" , only once, and it jumps back to the state "on" before 
being trapped again in the state "off" . This property ensures the renewal property. In conclusion, superstatistics is 
a form of very slow modulation, whereas renewal is compatible with a theory based on non-exponential distribution 
densities emerging from the superposition of many exponential waiting time distributions, provided that the system 
does not adopt the same exponential well for too many successive trapping events. 

As an example of physical process that might be conveniently described by the slow modulation perspective, we 
quote the interesting recent papers and [2(|. These authors fin d 1191 t hat their experimental result cannot be 
described by the two-state semi-Markov models proposed in Refs. [2l|, |22], being therefore incompatible with the 
renewal perspective. The chronological ordering of the off waiting times suggests the occurrence of a modulation 
corresponding to molecular conformational changes. 

It is therefore convenient to develop a technique of analysis of the experimental data that might help the investigators 
to establish the real nature of the process under study, namely whether renewal, slow modulation, or an intermediate 
condition, applies. This is the main purpose of this paper. We shall show that, although renewal and superstatistics 
can be used to produce the same waiting time distribution density, the former approach to complexity generates 
renewal aging, while the latter does not. The aging technique has already been adopted by Brokmann et al. |23| to 
prove that the physics of BQD is characterized by renewal. In this paper, in addition to confirming the conclusions 
of this earlier work, using the same aging technique, we make also an attempt at assessing if a condition intermediate 
between renewal and very slow modulation (superstatistics) might exist. 

The outline of the paper is as follows. In section[H]we show that renewal rests on resetting the system's memory after 
any event (collision). In Section [IHI we define modulation as a condition of time dependent rate, with no renewal. We 
explain the form of modulation adopted in this paper, and we show that in the fast condition it becomes compatible 
with non-Poisson renewal. We devote Section llVI to reviewing the phenomenon of renewal aging and we explain why 
the condition of very slow modulation should yield no aging. In Section with the help of artificial sequences, we 
illustrate the aging experiment that we propose for the analysis of experimental sequences. At the same time we 
explore the unknown region between renewal and very slow modulation. We examine and discuss real experimental 
data in Section IVll Finally, we devote Section IVTTl to concluding remarks. 



II. RENEWAL 



Let us consider a two-state renewal process, and, for simplicity, let us assume that the distribution of sojourn times 
in the state "on" is the same as the distribution in the state "off" . This assumption will prevent us from discussing 
with our techniques the interesting effect recently discussed by Verberk et al. [24), These authors discussed the case 
where both distributions have an inverse power law form with different indexes, (J-off and fi on and found that a Gibbs 
ensemble of trajectories moving from the beginning of the state "on", produce a fluorescence intensity decaying in 
time as a function proportional to l/£^°//~^°™, with > fi on . We leave the discussion of this interesting effect as 
a subject for future work. 

We assign to the Survival Probability (SP) of this process, ^(t), the inverse power law form 

with (j, > 1. This corresponds to the joint action of the time dependent rate [25| r(t) = r /(l+rit), with r = (p—l)/T 
and r\ = 1/T, and of a resetting prescription. To illustrate this condition, let us imagine the random drawing of a 
number from the interval / = [0, 1] at discrete times i = 0, 1, 2.... The interval I is divided into two parts, I\ and -?2j 
with Ii ranging from to pi, and I2 from pi to 1. Note that Pi = 1 — <ft < 1 and % << I, and, as a consequence, the 
number of times we keep drawing numbers from I\, without moving to I2, is very large. 
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Let us evaluate the distribution of these persistence times, and let us discuss under which conditions we get the SP 
of Eq. ft|. The SP function is the probability of remaining in I\ after n drawings, and is consequently given by 

n 

9{n) =l[ Pi . (2) 
i=i 

Using the condition q. L « 1, and evaluating the logarithm of both terms of Eq. @, we obtain: 

n 

M*W) = -E?i- (3) 

i=l 

The condition qi < < 1 implies that i and n of Eq. (|3J) are so large as to make qi virtually identical to a function of the 
continuous time t, q ?; = r(t) = r a r/(t), with rj(t) — 1/(1 + r\t). Therefore r(t) is a time dependent rate, with resetting 
as we hereby shall see. Thus, Eq. © yields the SP of Eq. QJ, and the corresponding waiting time distribution 
density, i/j(t), reads 

T-r/t-l 

W = 0»-^. (4) 

We denote as collisions the rare drawings of a number from li , followed by resetting. Thus the collisions occurring 
at times t\, Ti 4- t 2 , yield: ??(£) = 1/(1 + rii),0 < i < ri; 7y(i) = 1/(1 + T\{t — ti)),ti < t < n + t 2 , and so on. 
Note that 77(0) = 1 means that we prepare the system at time t = 0. We might adopt a coin tossing prescription to 
decide whether to keep or to change sign, after any collision. However, in this paper, as earlier pointed out, we do 
not pay attention to the problem of fluorescent intensity changing in time as an effect of ensemble average. Thus, for 
simplicity, our theoretical remarks refer to a sequence {r^}, where the times Tj are randomly drawn so as to yield the 
analytical form of Eq. (£1J, without assigning to them either an "on" or an "off" symbol. 



III. MODULATION 



The renewal condition described in Section ITT1 must not be confused with the case of a time dependent rate, with 
no renewal. The time dependence of r(t), with no renewal, might obey a deterministic or a stochastic prescription. 
An example of the former case is q(t) = A + Bcos(u)t), with no renewal. We think that the physical process studied 
by the authors of Refs. 19, 20], might be adequately described by a prescripion of this kind, not necessarily periodic, 
or quasi periodic, reflecting however the molecular conformational changes in time. 

The specific cases discussed in this paper are closer in spirit to the condition of stochastic dependence on time. 
This means that r(t) is a stochastic function of time, so that we have to study: 



=< exp (- J dt'r(t')\ > . (5) 

An interesting example of treatment of this kind is offered by the recent work of Brown j2|J . It has to be pointed 
out that the evaluation of the characteristic function of Eq. (J5J might be a difficult problem, but in the limiting cases 
of very fast or very slow modulations. The first condition departs from the non-Markov condition of interest for us in 
this paper. The latter condition does not, and can be adopted to derive the non-exponential behavior of ^(t). 

Let us assume that the fluctuation r(t) has an equilibrium distribution, p eq (r). In the special case where the time 
scale of the fluctuation r(t) is virtually infinite, the SP ^(t) becomes: 

/•oc 

*(t) = / drp(r)exp(—rt). (6) 
Jo 

In fact, let us consider a Gibbs ensemble of identical systems, obeying Eq. JSJ). At the moment when the observation 
begins, at time t = 0, each of these systems has a rate r, given by the distribution p eq (r). If the time scale of the rate 
fluctuations is virtually infinite, during the observation process each system will keep unchanged its own rate, thereby 
producing Eq. ©. This makes this picture essentially identical to the approach to complexity recently proposed by 
Beck H3. 
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According to a mathematical prescription borrowed from Beck [23 , we find that the analytical expression of Eq. 
is recovered by using the following form: 

Peg(r) = T ^_ l) r»- 2 exp{-rT). (7) 

This proposal was used in a later paper Q to develop a new approach to complexity, denoted as superstatistics. This 
approach is attracting the attention of an increasing number of researchers (see, for instance, |28l I29I H(j), and for 
this reason is worth of consideration. 

Note that once the inverse power law form of Eq. for the SP 'J'(i) is realized, the corresponding waiting time 
distribution density is given by 

poo 

"4>(t) = / drp(r)rexp(—rt). (8) 
Jo 

In principle, the adoption of the time-dependent rate prescription of Section |H] would make it possible to describe 
modulation processes of any kind. It would be enough to set r\ — and replace tq with q(t), and make it either 
deterministic or stochastic. However, to generate the time series {t{\ corresponding to the modulation prescription, 
fast, slow or of intermediate speed, we adopt a different procedure. We select from the distribution of Eq. J7J the 
sequence {rj}. For any exponential waiting time distribution ipj(t) = rjexp(-rjt), we select times, t(, ....,t^ . 

A more realistic picture might allow to fluctuate. We expect, however, that for ordinary fluctuations about a 
common value < N} n >= N m , the physics of the process does not change. Thus, for simplicity we assign to all 
the numbers the same value N m . The time series to compare to the one derived according to the prescription 
of Section [H] is defined by {t^} = t-J, ....t^ ,rf , ....t^ , .... The benefit of this criterion is that N m — 1 makes the 
resulting time series equivalent to that generated by the models of Refs.p^L IT6l UtI and, consequently, to the 
renewal prescription of Section UTI 



It is important to stress that the ideal condition of totally renewal process and the ideal condition of infinitely 
slow modulation are characterized by a marked difference concerning ergodicity. This important issue has been 
recently discussed by Margolin and Barkai [3l|]. In the totally renewal case, with /1 < 2, the system does not admit 
any stationary condition |32| and the stationary correlation function does not exist. The non-stationary correlation 
function can be defined making an ensemble average, as illustrated, for instance, by the authors of Ref. 33]. The 
authors of Ref. , on the contrary, adopt a single trajectory picture and study the time averaged intensity correlation 
function of the BQD signal, supposed to be totally renewal, as well as with identical "on" and "off" distribution, 
as assumed in this paper. The correlation function is a stochastic property characterized by U- and W-shaped 
distributions. 

The case of infinitely slow modulation would suggest the adoption of an ensemble rather than individual trajectory 
treatment. However, in this paper we adopt the individual trajectory treatment also for the case of slow, but not 
infinitely slow, modulation. It is expected that in this case ergodicity is not violated, in a striking contrast with the 
condition of total renewal [3l| . 

We think that moving from N m = 1, where the properties found by Margolin and Barkai |^| apply, to N m — 00, 
where only the ensemble treatment is possible, implies the exploration of an unknown region, of which this paper 
affords a preliminary treatment. For the reader to appreciate this aspect, we would like to introduce the concept of 
pseudo event. This concept is similar to that proposed in an earlier publication Q. The authors of this paper [34| 
found that in some problems of medical interest the connection between scaling and waiting time distribution does 
not correspond to the prescription of the renewal theory This is so as a consequence of the fact that the times of the 
series under study turned out to be correlated [34| . 

In the case of modulation, we define as pseudo events all the drawings of waiting times from the same Poisson 
distribution, after the first drawing. In the case N m = 1 there are no pseudo events. In the case N m — 2 there 
is one pseudo event, and so on. The quantity N m — 1 defines the number of pseudo events per critical event. By 
critical event, we denote the drawing of a given Poisson parameter r. In practice, the occurrence of a critical event 
corresponds to the first drawing of a waiting time from a Poisson distribution with rate Tj , namely the time t[ . 

The drawing of the next waiting times from the same distribution, implies a subtle deviation from renewal. This 
form of correlation is not easy to detect. In fact, although consecutive sojourn times are drawn from the same Poisson 
distribution, they arc by definition independent one from the other. If the time correlation function is < TjTj >, with 
f = r— < t > for a finite portion of the sequence, we expect it to yield (< t 2 > — < r > 2 )Sij, with Sij is the delta 
of Kronecker. 

It is a striking and surprising fact that the correlation produced by modulation is invisible to the ordinary correlation 
test. This is so because the sojourn times, although derived for an extended period of time from the same Poisson 
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distribution, are randomly drawn. The aging experiment reveals this subtle form of correlation. We refer to the times 
that are correlated the ones to the other as pseudo events, regardless the origin of this correlation, which might occur 
in the form discussed in the earlier publication [^H or in the even more subtle form of this paper. We think that the 
aging reduction depends on the ratio of pseudo to critical events, regardless the origin of correlation. In this paper we 
find that the aging reduction depends indeed on the ratio of pseudo to critical events, of the type here introduced. On 
the basis of this result with artificial sequence, we make also an attempt at evaluating the amount of pseudo events 
that might be present in the real BQD time series. 



IV. AGING AND MODULATION 

For an intuitive description of the concept of aging we shortly review the treatment of the earlier work of Ref. |35| . 
The renewal condition of Section [n] can also be realized with the following dynamical model. A particle moves in the 
interval / = [0, 1] driven by the equation of motion 

§-»•■• (9) 

with < a << 1 and z > 1. When it reaches the point y = 1 it is injected back in a random position between and 1 
with uniform probability, thereby producing another extended time of sojourn within the interval /. The connection 
between the waiting time distribution density and the uniform initial distribution is given by 

ip(t)dt = p(y )dy . (10) 
It is easy to prove that the resulting time distribution density is given again by Eq. (0J with 

z 

H = 7 (11 

z — 1 



and 



T=^— i. (12) 



The uniform back injection is equivalent to the resetting prescription of Section ^ and, in fact, this model is 
renewal, and it is equivalent to the model of Section^ but its adoption in this section serves better the purpose of 
explaining renewal aging and the lack of aging in the case of very slow modulation. 

The waiting time distribution density given by Eq. l|lUfl corresponds to beginning the observation at the moment 
when the system is prepared in the uniform distribution p(yo) = 1. As a result of the injection back process this 
distribution changes upon time change. If the observation of the first times of sojourn is made at a later time t a > 0, 
the corresponding waiting time distribution density is given by 

ip ta (t)dt = p(y , t a )dy . (13) 

The dependence of tpt a (t) on t a is the renewal aging that we want to assess in this paper by means of a suitable 
numerical experiment. An exact expression for V'taM is available |35l l36||. but, since it is not expressed as a simple 
analytical formula, it is not suitable for the practical purposes of this paper. For this reason, we prefer to adopt the 
expression: 

Mt)= & dV t {t + V \ (14) 

where K± a is a suitable normalization constant. The meaning of this approximated expression is evident. We assume 
that the first sojourn times observed might have begun prior to t — t a , anywhere between t = and t = t a , with the 
restrictive condition that the earlier laminar region, only one, began at t = 0. Actually, the last laminar region might 
be at the end of a sequel of an arbitrarily large number of jumps, thereby generating corrections to Eq. I|14|) . In Ref. 
|37| the accuracy of this approximation, in the case of inverse power law waiting time distributions, was examined, 
and found to be very good. In this paper we shall make a discussion of the key results on the aging experiment on 
BQD systems, taking into account the error associated to this approximated formula. 

It is possible to predict that the case of very slow modulation does not yield aging. In the case of a modulated 
Poisson process, we replace the model of Eq. © with 



(15) 
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This means that we set 2=1 and we replace the parameter a with the time dependent rate r(t). The equation of 
motion for p(y, t) is given by 



d_ 

Of 



p(y,t) = r(t) 



d_ 

dy 



yp(y,t) +p(±,t) 



(16) 



Note that the second term of on the right hand side of this equation corresponds to the back injection of the particle, 
when it reaches the border y = 1, and thus to the resetting process of the renewal model of Eq. ©. When r(t) does 
not depend on time, Eq. I|16|l represents a Poisson process. Let us focus our attention on the case where r(t) is a 
stochastic function of time. If it is very fast, r(t) must be replaced by < r{t) >, the process becomes Poisson again 
and it departs from the modulation model adopted in this paper (see Section ITTTI . Anyway, in accordance to a well 
known notion, this Poisson process does not yield aging. The model of Eq. Ijl6(l becomes equivalent to the modulation 
model of this paper when the fluctuation of r(t) is very slow. We see that the equilibrium distribution coincides with 
the initial flat distribution. Thus, we cannot adopt the departure from the initial distribution as a way to define the 
system's age. In the case of a virtually infinitely slow modulation, the system leaves for a virtually infinite time in a 
Poisson condition, with no aging whatsoever. 



V. AGING EXPERIMENT ON ARTIFICIAL DATA 



This section is devoted to illustrating, with the help of artificial data, a technique of analysis aiming at a quantitative 
evaluation of the degree of renewal properties of a given time series. We refer to this kind of analysis as aging 
experiment. 

It is important to notice that the analysis of real data implies the observation of only one single sequence. In this 
case we must turn a single sequence into a very large number of sequences of the same age. The first sequence is 
the sequence, artificial or experimental, to analyze, beginning at time t = with the system being located at the 
beginning of a state, either "on" or "off". The second sequence is obtained from the first, canceling the first state, 
namely, shifting the first sequence towards the time origin by the quantity equal to the time duration of the first 
state, so that the second sequence begins at time t = 0, when the system begins sojourning in the second state of 
the first, or original, sequence. On the same token, the third sequence begins at time t = 0, when the system begins 
sojourning in the third state, and so on. Thus the waiting time distribution ip(t), t a = 0, is the distribution of the 
time durations of the first states. To do the aging experiment we set t a > and we record the time lengths of the first 
states observed in that time position. With this prescription we define ip ta {t)- A quantitative definition of amount of 
aging is more properly done using the SP, defined by 

/oo 
iptAt'W, (17) 

rather than ifj ta (t). 

From now on we denote by ip exp (t) and if? exp (t) the waiting time distribution densities and the SPs, respectively, 
derived from the experimental data. Of course, in the case of artificial data, where the sequence is realized for the 
specific purpose of producing the function 'i'(t) of Eq. 0) and ip(t) of Eq. (@J, the experimental functions coincide 
with the corresponding theoretical prescriptions. Then, we define the corresponding aged distributions using Eq. I|14f) . 
This expression is not exact. However, it is convenient for the purposes of this paper, where the experimental error 
is expected to be larger than the discrepancy between Eq. I|14fl an d the exact prescription. Then, we denote with 
\&j en (i) the SP derived from the experimental observation, namely from *f?Q Xp (t) = fyl x = (t), by means of Eq. I|14fl . 

The prescription adopted in Section ITTTI to produce the time series {t^} with a changing N m for N m — > oo becomes 
coincident with the slow modulation of Eq. (|16l) . Thus, we expect no aging in this limiting case. In the opposite limit 
with N m — 1, the sequence is renewal. Thus, we expect the maximum amount of aging. In other words, the renewal 
condition should yield 

K p (t) = *IT(*), ( 18 ) 

whereas the condition of very slow modulation should produce no aging, a property described by: 

= % en (t) = % xv . (19) 

Eqs. I|18fl an d (|19fl refer to two limiting conditions: the condition of Eq. i|18|) corresponds to total aging, with no 
pseudo events, while the condition of Eq. I|19(l stems from a process dominated by Poisson pseudo events, with a total 
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FIG. 1: The SP *^ p (t) as a function of time. The three figures Fig. [T(a)1 [T(b)1 and [T(c)| refer to N„ 
N m — 100, namely to 0, 9 and 99 pseudo events per critical event, respectively. The curves of each figure are obtained from 
artificial sequences yielding for the waiting time distribution density the same inverse power law form of Eq. with jj, — 1.8. 
For each figure the three distinct ages t a = 0, 20, 60 (from the bottom to the top) are considered. The renewal predictions and 
the modulation results are denoted by full lines and circles, respectively. 



lack of aging. It is important to point out that the results presented in Section [VII do not involve any assumption on 
the form of the waiting time distribution. 

Let us introduce here another important ingredient of our analysis, the aging intensity function: 



w a n (r)-n xp (ry 



(20) 



Eqs. (|18fl and 119fl yield J a (r) = 1, and I a (T~) = 0, respectively, thereby indicating that Eq. (|2U[1 is a proper aging 
intensity indicator, with 1 and representing total aging and lack of aging, respectively. In principle, this function 
should decrease from 1 to upon increase of the number of pseudo events. 
The numerical results of Figs. 1(a) 



to 



1 re- 



confirm the expectation that the larger the number of pseudo events, the 
smaller the aging intensity. In fact, in accordance with the earlier theoretical remarks, we see that in the case of no 
pseudo event, N m — 1, reported in Fig. |l(a)j modulation and renewal yield the same amount of aging. The occurrence 
of 9 pseudo events, n amely the case N m = 10 illustrated in Fig. 1(b) is already enough to significantly reduce 
modulation aging. Fig. |l(c)| shows that N m = 100 yields an even larger aging intensity reduction. In conclusion, these 
numerical results confirm the theoretical expectation that the infinitely slow modulation, N m = oo, should produce 
no aging. 

The adoption of the aging intensity function of Eq. I|20|) allows us to express the aging reduction illustrated by 
the earlier figures in a quantitative way. The analysis of the numerical simulations of both renewal and modulation 
models shows that the aging intensity function I a (j) for t — > oo tends to an asymptotic value, / a (oo). This property 



is shared by the real data analyzed in Section rVTl We estimate this value and we use it, in both the case of artificial 
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FIG. 2: The aging intensity indicator I a (oo) as a function N m , namely iV m — 1 pseudo events per critical event, at different 
age, namely different values of t a . The aging intensity is defined by Eq. (1201 . These results refer to artificial sequences with 
fj, = 1.8. Physical conditions with different ages and the same number of pseudo events collapse into the same 7 a (oo), thereby 
becoming indistinguishable in the scale of this figure. 



sequences of this section and in the case of real data of Section IVT1 to define the time asymptotic intensity of the aging 
indicator. 

Fig. illustrates the application of this procedure to the case of artificial sequences. We see that the aging intensity 
decreases with the increase of the number of pseudo events per critical event. We fit I a (oo), as a function of N m , with 
the inverse power law (N m )~ a , with a = 0.70 ± 0.02. The aging intensity indicator, 7 (co), which should hold the 
value of 1 when there are no pseudo events, actually slightly exceeds this value with no pseudo event and decreases 
by a factor of 10, with increasing the number of pseudo events per critical event from to 99. In Section IVTI we shall 
refer to Fig. [21 to estimate the amount of pseudo events per critical event present in the real BQD data. 

The aging intensity overestimation, with no pseudo events, is a consequence of the fact that Eq. (|14|l is not exact. 
Let us study the effects of this inaccuracy by means of artificial sequences. We use artificial sequences derived from 
the renewal prescription, namely, by random drawings of numbers from the inverse power law distribution of Eq. 0} , 
yielding a SP of the form of the Eq. (JJJ. We use the values fj, = 1.65 and \i = 1.8, which are typical values of the "off' 
sequences studied in Section IVll We make a Monte Carlo simulation and we produce curves of the kind illustrated 
in Fig. |3| showing the time dependence of the aging intensity function. We see that I a (T) becomes virtually time 
independent after a short transient. Let us notice that the levels of these plateaus exceed the maximum value of 1. 
However, we see that this level becomes closer to 1 with the age increase from t a = 20 to t a = 220. Of course, the 

1.4 
1.35 
1.3 
1.25 
(I«)(t) 1.2 
1.15 
1.1 
1.05 
1 

500 1000 1500 2000 2500 3000 



FIG. 3: The aging intensity J Q (r) of Eq. I20H as a function of r. We study the case of renewal artificial sequences with an 
inverse power law form, with /i = 1.8. We study the time evolution of this indicator at two different ages, t a = 20 and t a = 220. 
The attenuation of the overestimation effect with the increasing t a is evident. 
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V = 


1.65 


M = 


1.8 


t a 


/ O (oo) 


&ren (-^a ) 


ia(oo) 


^ren( la ) 


20 


1.3734 


0.0011 


1.3378 


0.0047 


60 


1.31828 


0.00070 


1.2765 


0.0026 


100 


1.3007 


0.0012 


1.2608 


0.0017 


140 


1.29325 


0.00040 


1.2502 


0.0014 


180 


1.28764 


0.00036 


1.2420 


0.0014 


220 


1.28243 


0.00037 


1.2344 


0.0012 



TABLE I: This table summarizes the results of Monte Carlo simulations done to evaluate the aging intensity I a (oo). These 
results refer to the case of renewal artificial sequences, with an inverse power law form. The values of p adopted, /i — 1.65 and 
(i = 1.8, are the typical values of the real BQD sequences of "off' states. oy en is the standard deviation. 

property of a fast attainment of a plateau is shared also by the cases, not shown here, of a non vanishing number 
of pseudo events, with the plateau height significantly smaller than 1 for N m of the order of 10: the results of Fig. 
@ have been obtained by evaluating the level of this plateau. Applying the same criterion to the case of no pseudo 
events, we get the asymptotic values I a (oo), illustrated by Table [I] 

In conclusion, the aging indicator I a (oo) should be independent of the system age and equal to 1: in fact, we are 
analyzing artificial sequences corresponding to the renewal model. We see, on the contrary, that the aging intensity 
indicator exceeds the maximum value of 1. This is due to the inaccuracy of the formula of Eq. J20J), which is known 
|35| to become accurate for t a — > oo. In accordance with this fact, we see from Table [l]that the aging intensity 
overestimation tends to vanish with increasing values of t a - 

VI. AGING EXPERIMENT ON REAL DATA 

We are now in a position to analyze the experimental results on real BQD data sequences. The experimental 
data discussed in this Section have been obtained by Prof. M. Kuno and V. Protasenko, Dept. of Chemistry and 
Biochemistry, University of Notre Dame. Our data set consists of 32 sequences of BQD fluorescence intensities. Each 
sequence contains 1 hour of records, sampled every ms, for a total amount of 360000 data per sequence. A sample of 
the data studied in this section is shown in Fig. 




FIG. 4: A sample of the BQD data examined in this section. 



In order to separate the "on" from the "off" state it is necessary to define a threshold intensity, above which the 
signal corresponds to a " on" state and below which it signals the "off" state. We establish this separation along the 
lines adopted in an earlier work (see p^L namely with an iterative procedure for the search of a repartition that 
make it possible for the fluctuations of the two states not to intersect with the separation line. At the end of this 
iteration process the variances of the two states get a well defined value, with a denoting the variance of the "off" 
state: the threshold turns out to be located at the value 2a over the "off" state. 

After defining an alternate sequence of "on" and "off" states, we make the aging experiment with the criterion 
described in Seclion lTvl Actually, we make three different kinds of aging experiment. The first and second, considering 
the waiting times only of the "on" and "off" states, respectively. This means that we sew the beginning of "on" ("off") 
state to the end of the immediately preceding "on" ( "off" ) state. The third experiment is done on the whole sequence 
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of waiting times, with the jump from one state to the other signaling the presence of an event, whose statistical 
properties are studied regardless of whether it corresponds to jumping from the "on" to the "off" state, or vice-versa. 

A. Aging of the "on" state 

Here we discuss the results of the first kind of aging experiments, on the "light on" state. Fig. [S] shows an example 
of these aging experiments. This figure refers to a case where the aging experiment on the BQD sequence is done for 
several values of t a , ranging from 20 to 220 in steps of 40. The thin continuous line represents 1 $q XP , the dotted line 
refers to ^ xp , while the thick continuous line represents ^\ en . We see that in this case the condition of Eq. (f 1 81) is 
fulfilled with a very good accuracy. The systems ages and it does according to the non-Poisson renewal theory. In 
fact, aging implies that the system obeys a non-Poisson statistics, and the fulfillment of Eq. (|18fl means that this 
non-Poisson statistics is renewal. It is also worth pointing out again that this observation does not involve that the 
deviation from Poisson statistics is realized through inverse power laws, as assumed for simplicity in Section [H] In 
fact, the experimental waiting time distribution and SP are not inverse power law, or, at least do not correspond to 
an inverse power law with a well defined index. However, they depart from the exponential condition enough as to 
generate the aging effects illustrated by Fig. [5j 



Aging experiment for the "light on" state 



9- 




0.1 1 10 100 1000 

x 

FIG. 5: The SP ^(t) as a function of time. The experiment is carried using only the "light on" waiting time distribution. The 
thin continuous line indicate ^Q Xp (t), the dotted line refers to ^ exp (t), and the thick continuous line is f r en- From the bottom 
to the top the curves refer to t a with the values 20, 120 and 220. 



B. Aging of the "off" state 

Let us discuss aging for the sequence of waiting times in "off" state. A sample of these results is shown by Fig. 
Also in this case, the system is aging and its behavior is very well described by the renewal theory. 

By visual inspection, we can stress some differences between the result of the latter and the former experiment. 
The main difference seems to be that the distribution of the time of sojourn in "light on" state is truncated after 
about two decades, while the inverse power law distribution of the "off" states holds longer. This result confirms the 
earlier observation of Chung and Bawendi 5]. For this reason, it turns out to be difficult for us to estimate the index 
of the waiting time distribution of the "on" states. In the case of the "off" states, instead, we can do that, since the 
power law behavior is more distinct that in the earlier case. Our estimations of the power law exponent, in this case, 
ranges from 1.65 ± 0.02 to 1.80 ± 0.05. We note that the artificial sequences of Table I refer to these values. 
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Aging experiment for the "light off" state 




0.01 1 ' ' ' ' 1 

0.1 1 10 100 1000 

T 

FIG. 6: The SP as a function of time. The experiment is carried using only the "off" waiting time distribution. The thin 
continuous line indicate ^fg Xp (t), the dotted line refers to ^ exp (t), and the thick continuous line is f r en. From the bottom to 
the top the curves refer to t a with the values 20, 120 and 220. 

C. Aging of the "on-off" state 

Let us discuss now the result of the third aging experiment, done on the whole sequence, with the transitions from 
the "on" to the "off' states as the markers of the significant events to analyze with the aging experiment. Fig. [7| 
shows a sample of this third kind of aging experiment. Also in this case it is evident how the theoretical predictions 
of the renewal theory fit very well the experimental results. 

In order to support quantitatively the conjecture made on the basis of visual inspection of the results, in the next 
subsection we shall adopt the aging intensity indicator introduced in Sec. [V] 

D. Aging intensity 

We recall that the adoption of Eq. (|14fl yields for the aging indicator I a (oo) values larger than 1, namely, we find 
that Eq. H14fl overestimates the aging intensity (see Table QJ. With Fig. 2 we also found that a number of pseudo 
events of the order of ten significantly reduces the aging intensity. 

We can now use these indications for a rough estimation of the number of possible pseudo events present in the 
sequences of sojourn times in the "light on" state. To make the evaluation of the aging intensity value as statistically 
accurate as possible, we adopt a Gibbs ensemble average. In other words, we make an average over all the sequences 
at our disposal, after assigning to them the same age t a . This is done by preparing all the sequences in such a way 
that at t = each of them begin at the beginning of the time of sojourn in the "on" state, the "off" state, or the 
"on" or "off" state, according to the kind of experiment under study, of the first, second and third type, respectively. 
Then we set for all sequences the beginning of the observation process at the same time t a . We obtain the mean aging 
indicator, indicated by < I a >. We do the same experiment for different values of t a - Fig.|S]shows the time evolution 
of I a (i~) for a given sequence of "light on" states corresponding to that used to derive the results of Fig. [31 The aging 
indicator I a (oo) for any sequence of this kind is obtained by making an average on r, and the fluctuations around 
this mean value are used to define the measurement error. 

The average on all the sequences of the Gibbs ensemble are used to obtain the values of < I a > reported in the 
tables. Tables ITT1 IIHI and IIVI report the value of < I a > for the "on", "off" and "on-off" waiting time distributions, 
respectively. In these tables, a represents the maximum likelihood estimation of the standard deviation of < I a >■ 

These results show that the aging of both the "off" and "on-off" waiting time distributions are very well described 
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Aging experiment for the "on-off" state 




0.01 1 ' ■ ■ ' 1 

0.1 1 10 100 1000 

T 

FIG. 7: The SP as a function of time. The experiment is carried using both the "on" and "off' waiting time distribution. 
The thin continuous line indicate ^g Xp (t), the dotted line refers to 0! exp (t), and the thick continuous line is 'I'ren. From the 
bottom to the top the curves refer to t a with the values 20, 120 and 220. 




2 50 5DO Vf50 1O0O 1^SO 1 750 POOO ^Pf50 

FIG. 8: Example of the aging intensity function I a (r) for the "light on" state. 



ta < la > a {I a) 

20 0.676 0.026 

60 0.690 0.023 

100 0.810 0.021 

140 0.795 0.020 

180 0.740 0.022 

220 0.868 0.019 



TABLE II: This table shows the intensity of aging Eq. 1201 for the "light on" state, for different values of t a , averaged over the 
ensemble. 
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ta 


< la > 


°%) 


20 


1.040 


0.020 


60 


0.992 


0.012 


100 


1.030 


0.014 


140 


1.010 


0.014 


180 


0.947 


0.013 


220 


0.962 


0.015 



TABLE III: This table shows the mean aging intensity of Eq. (1 li ( ) tt for the "off" state, at different values of t a . The average is 
carried out on the whole experimental sequences. 



ta 


< la > 




20 


0.949 


0.013 


60 


0.894 


0.011 


100 


0.914 


0.011 


140 


0.865 


0.011 


180 


0.823 


0.012 


220 


0.810 


0.011 



TABLE IV: This table shows the mean aging intensity Eq. 1201 for the "on-off" state, at several values of t a . The average is 
carried out on all the experimental sequences. 

by the means of the renewal theory. The accuracy of the renewal prediction, become worse for the "light on" state, 
especially for low values of the parameter t a , and it improves with the increase of t a . By comparing these values to 
the curve of Fig. 2, and taking into account that with N m = 1 the aging intensity indicator overestimates the aging 
intensity of the artificial sequences, we cannot rule out the possibility that the "light on" state might involve the 
presence of from 5 to 6 pseudo events. 



VII. CONCLUDING REMARKS 

The main result of this paper is the proof that BQD data obey non-Poisson renewal with a good accuracy. This 
complexity condition does not stem from modulation (superstatistics). The statistics of the "on" states is not identical 
to the statistics of "off" states. The latter case is closer than the former to the inverse power law picture of Section 
IITI However, both waiting time distributions depart significantly from the exponential form and produce significant 
aging effects. 

We confirm the general opinion that the BQD phenomenon obeys renewal prescription, leaving open, however, the 
possibility that a finite amount of pseudo events might be involved, especially for the "light on" state. This suggestion 
emerges from the results of Tables Inland IIIII compared to those of Fig. 

We have to point out that the field of sin gle- system spectroscopy is wide and there are examples of process where 
slow modulation conditions seem to apply 0, |2(], HU . This paper affords prescriptions of statistical analysis that 
might turn out to be useful to study the unknown territory between totally renewal and infinitely slow modulation. 
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